Topological insulators in filled skutterudites 
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We propose new topological insulators in cerium filled skutterudite (FS) compounds based on ab 
initio calculations. We find that two compounds CeOs4Asi2 and CeOs4Sbi2 are zero gap materials 
with band inversion between Os-d and Ce-f orbitals, which are thus parent compounds of two and 
three-dimensional topological insulators just like bulk HgTe. At low temperature, both compounds 
become topological Kondo insulators, which are Kondo insulators in the bulk, but have robust 
Dirac surface states on the boundary. This new family of topological insulators has two advantages 
compared to previous ones. First, they can have good proximity effect with other superconducting 
FS compounds to realize Majarona fermions. Second, the antiferromagnetism of CeOs4Sbi2 at low 
temperature provides a way to realize the massive Dirac fermion with novel topological phenomena. 

PACS numbers: 71.20.-b,73.43.-f,73.20.-r,75.30.Mb 
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The topological insulator (TI) is a new state of quan- 
tum matter attracting extensive attention for both 
fundamental physics and technical applications recently. 
TIs have insulating energy gap in the bulk and accommo- 
date gapless edge or surface states which are protected 
by time-reversal symmetry (TRS). With TRS broken on 
the surface of a three-dimensional (3D) insulator, the ef- 
fective electromagnetic response is described by the topo- 
logical term S = (9/2w)(a/2n) J d 3 xdtE • B, with a the 
fine-structure constant and the parameter 9 = or ir for 
trivial insulator and TI, respectively |5|. This topolog- 
ical response supports many novel topological phenom- 
ena, such as image magnetic monopole induced by a point 
charge @, topological Faraday and Kerr effects [Bf and the 
realization of axion field in condense matter physics 0. 
Therefore magnetically doped TI for magnetic impurities 
and ferromagnetism effects are of particular interest. In 
addition, superconducting proximity effects on the sur- 
face states of a 3D TlQ have been proposed to realize 
the Majorana fermion which is a promising candidate 
for quantum computation application. In this letter we 
report new TI materials which could realize both the an- 
tiferromagnetic (AFM) TI and the Majorana fermion. 
Furthermore they are found to be topological Kondo in- 
sulators at low temperature. 

HgTe quantum wells were the first TI pre- 
dicted theoretically and subsequently observed 
experimentally [l^- The basic mechanism of band 
inversion driven by spin-orbit coupling (SOC) was 
discovered in this work and provides a template for most 
TIs discovered later. Thereafter, many TIs have been 
proposed and experimentally measured, e.g. Bi2Se3 and 
Bi 2 Te 3 crystals and TlBiTe 2 0, (H crystals. 
Most of them are known as thermoelectric materials 
before. At the same time theoretical calculations have 
predicted many other TI materials which are waiting 
for experimental verification 0, Here we find that 
CeOs4Asi2 and CeOs4Sbi2 within the filled skutteru- 
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FIG. 1: (Color online) (a) Crystal structures in bec lattice for 
filled skutterudites RT4X12 {R— rare-earth, T= Fe, Ru or Os, 
and X— P, As or Sb) in a bec lattice and (b) corresponding 
first Brillouin zone, in which F, N and H are time-reversal- 
invariant points. 



dite(FS) class are TI materials. The FS compounds [16| 
have a chemical formula RT4X12 (R= rare-earth, T= 
Fe, Ru or Os, and X— P, As or Sb), in which heavy 
elements are expected to induce strong SOC. Similar 
to Bi2Se3 materials, they are also known for excellent 
thermoelectric properties [T3, IH| Moreover, they exhibit 
a rich variety of electronic and magnetic ground states 
at low temperature [13- 21 1, including superconductivity, 
ferromagnetism, antiferromagnetism and Kondo insu- 
lator behavior with hybridization gaps. Among these 
compounds, most of Ce-based FSs are reported to be 
insulators [Til. All these reasons motivate us to 

investigate the Ce-based FSs as TI candidates which 
can accommodate magnetism and superconductivity 
proximity. 

The FSs crystallize in the bec structure (Im3) shown 
in Fig.l and have inversion symmetry. The R atoms, 
the inversion centres, are located at the body centre 
and corners of the cubic structure and surrounded by 
a cage of corner-sharing T4X12 octahedra. We calcu- 
lated the electronic properties of CeTiX\2 by density- 
functional theory (DFT) within Perdew-burke-Ernzerbof 
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TABLE I: The calculated band gap and the Z2 topological 
index. We listed the calculated band gap _E (calc.) together 
with the experimental value -E 9 (exp.)(see refl2lPl in unit of eV. 
The products of parities for time-reversal invariant points in 
the bcc Brillouin zone are shown, including the F point and 
one H point and six A" points. 



FSs 


i? 9 (calc.) 


E g (cxp.) 


Parity 


Z 2 


CeFe 4 Pi2 


0.38 


0.13 






CeFe4Asi2 


0.16 


0,0.01 






CeFe 4 Sbi2 


0.08 





r(-) 




CeRu 4 Pi2 


0.12 


0.086 


H(-) 


(0;000) 


CeRu4Asi2 


0.13 


0,0.0043 


iV(-) 




CeRu4Sbi2 


0.12 









CeOs 4 Pi2 


0.12 


0.034 






CeOs4Asi2 


0.00 


0.0047 


r(+) 

H(-) 


(1;000) 


CeOs4Sbi2 


0.00 


0.0009 


AT(-) 





type generalized gradient approximation(GGA).The Vi- 
enna ab initio simulation package (vasp) 23] with the 
projected augmented wave method are employed. We 
adopted the experimental lattice constants (see ref. flit ) 
and fully optimized the atomic positions in the unit cell. 
SOC is included in our calculations. In addition, wien2k 
package [iij] was also used for cross-checking. 

The calculated GGA band gap as well as the exper- 
imental values are shown in Table I. Our results are 
well consistent with previous calculations for CeFe4Pi2 
and CeFe 4 Sbi 2 [25l and CeOs 4 Sbi 2 [26j]. Most materials 
have a small gap while CeOs4Asi2 and CeOs4Sbi2 have 
a zero gap. We show the band structure of CeOs4Xi2 
in Fig. 2. Take CeOs4Pi2 as an example. The Ce-4/ 
states hybridize with both Os-5<i and P-3p states near the 
Fermi energy. A band gap opens due to their hybridiza- 
tion. The hybridization can be overestimated, since GGA 
tends to underestimate the inter-atomic correlation of 
Ce-4/ electrons. Then the calculated band gap is ex- 
pected to be overestimated [2^], which is different from 
the usual gap underestimation of DFT calculation. This 
can explain why our calculated band gap is larger than 
the experimental one. For all these three FSs in Fig. 2, 
the Ce-f states exists near the valence top and becomes 
dominant as the narrow spin-orbit splitting bands at the 
conduction bottom. 

In order to determine their topological features, we use 
the parity criteria proposed by Fu and Kane [53] to calcu- 
late the Z2 topological index. The Z2 index is determined 
by the parity of occupied bands on each time-reversal in- 
variant momenta. The bcc Brillouin zone has eight time- 
reversal-invariant points, including the T point, six N 
points equivalent to (7r,0,0) or (7r,7r,0) by point group 
symmetry, and one H point (7r,7r,7r). We listed the prod- 
uct of parities for all occupied states at these /c-points 
and corresponding Z2 index in Table I. Among them FSs 
CeOs4Asi2 and CeOs4Sbi2 are found to be topological 



nontrivial with Z2(l;000). Compared to those topologi- 
cal trivial FSs, e.g. CeOs4Pi2, a band inversion between 
Os-d (T^) and Os-f (Y^) bands is found in the band 
structures, as shown in Fig. 2. The d(r^") state enters 
the conduction bands. Both the valence band maximum 
and conduction band minimum at Y are / (Yq 7 ) states 
and touch with each other. The degeneracy induces a 
zero band gap. This is similar with HgTe which has s 
{Te)-p (Tg) band inversion. In order to understand the 
band structure better, we showed in Fig. 3 the Fermi sur- 
face at Ep—0.0 and -0.1 eV cases for CeOs 4 Sbi 2 . At the 
Fermi surface, there is a small electron pocket around H 
point from / bands and a tiny hole pocket at T. For Ep 
slightly below zero, there is a hole packet around Y with 
a little anisotropy, which is due to the top valence band. 

To drive CeOs 4 Asi2 and CeOs4Sbi 2 into real topolog- 
ical insulators, the degeneracy at Y point needs to be 
lifted. There are several possible ways to lift the de- 
generacy. Firstly, it is possible to fabricate quantum 
wells using CeOs4Asi 2 or CeOs4Sbi2 together with above 
topological trivial FSs as barrier layers, like HgTe-CdTe 
quantum wells [tl [l(|. Then the 2D TI (quantum spin 
hall effect) can be realized. Secondly, strain can lift up 
the r^" 7 degeneracy and open a a gap at Y. Then these 
two materials will become 3D topological insulators, sim- 
ilar to strained HgTe 27-29]. For example, 2% uniaxial 
strain along (001) direction can open an energy gap at 
r for CeOs 4 Asi2 in our calculation, as shown in Fig. 
2(d). Thirdly, a new feature in these two materials com- 
pared to previously proposed TI materials is that they 
can be 3D topological Kondo insulators [13 ■ At low tem- 
perature (T < 135 K for CeOs 4 Asi 2 and T < 50 K for 
CeOs4Sbi2_), both materials are reported to be Kondo 
insulators |3li |32| . Though the GGA band gap is zero, at 
low temperature the residual carriers in the system can 
form Kondo singlets with localized / states, leading to a 
Kondo insulator behavior in transport. Moreover, both 
materials have p-type carrier 33J, l34| , so that the Fermi 
level does not cross the bands involved in the band inver- 
sion. Consequently, when the carriers become insulating 
due to Kondo effect, the topological surface state due to 
the band inversion are expected to remain robust on the 
surface, so that the system becomes a "topological Kondo 
insulator" rather than ordinary Kondo insulators. On 
the contrary, if the carrier is n-type, the bands involved 
in band inversion will be buried below the Fermi level, 
so that the topological nontrivial nature of the system is 
lost for high enough carrier density. If the carrier density 
in this system can be tuned, we expect to see a topo- 
logical phase transition between the topological Kondo 
insulator phase for p doping and the ordinary Kondo in- 
sulator phase for n doping. It should be noticed that the 
topological Kondo insulators studied here are different 
from those proposed by Dzero et al[3(|, since the resid- 
ual carrier density in the FS compounds is low, so that 
the Kondo effect leading to a Kondo insulator state is 
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FIG. 2: (Color online) Band structures of (a) CeOs4Pi2, (b)CeOs4Asi2, (c) CeOs4Sbi2 and (d) strained CeOs4Sbi2. Red 
dots stand for the components of Ce-f states, black dots for Os-d states and green for p states. The size of dots represent the 
relative amplitude of corresponding components. The Fermi energy is set to zero. 



expected to be induced by magnetic impurities, rather 
than a Kondo lattice formed by /-electrons. On a first 
view the Ce filled FSs have similar electronic structure 
like the rare earth containing Heusler compounds, which 
are also predicted to be HgTe-like TIs recently (35l. [36j| . 
Among them the heavy-fermion behavior (e.g. YbPtBi) 
and superconductivity (e.g. LaBiPt) are closely related 
to / electrons. However, the topological feature of the 
Heusler compounds are related to a s — p band inversion, 
different from these FSs. 

Recent experiments reveal that CeOs4Sbi2 has a AFM 
phase below 1 k[37]]. (It should be noted in above calcu- 
lations we constrained the magnetic moments as zero to 
keep the TRS, since the total energy difference between 
with and without magnetic moments is negligible small 
and beyond the accuracy of DF calculations) The exis- 
tence of AFM order provides a possible way to realize the 
massive Dirac fermion state with many exotic topologi- 
cal phenomena without requiring magnetic doping. For 
example, the spin wave in the AFM phase is expected to 
have "axionic" coupling to photons, leading to the ax- 
ion polariton effect and a new kind of tunable optical 
modulatorQ. Moreover, in proximity with their super- 
conducting FS neighbors, e.g. \j&T/±X\2, both CeOs4Asi2 
and CeOs4Sbi2 can help to realize the topological super- 
conductivity as both 3D and 2D TIs. 

In summary, we proposed two FS compounds 



CeOs4Asi2 and CeOs4Sbi 2 as new topological insula- 
tors. Up to now the experimentally measured TI ma- 
terials include only s-p band inversion such as HgTc 
and p-p band inversion such as Bi2Se3 family. Differ- 
ent from them, these two TIs have d-f band inversion. 
They have two advantages compared to previous ones. 
Firstly, they can have good proximity with other super- 
conducting FS compounds to realize Majarona fermions. 
Secondly, the antiferromagnetism of CeOs4Sbi2 at low 
temperature provides a way to realize the massive Dirac 
fermion with exotic topological phenomena. Moreover, 
their rich physics offers a platform for topological Kondo 
effects. At low temperature, both compounds become 
topological Kondo insulators, i.e. they are Kondo insula- 
tors in the bulk, but have robust Dirac surface states on 
the boundary. After the completion of this work, we be- 
came aware of recent work on the skutterudite CoSb3 [IH , 
which is near the topological critical point. 
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(a)E F =0.0 (b)E=-0.1 



FIG. 3: (Color online) Fermi surface for CeOs4Sbi2 at Ef = 
(a) 0.0 and (b) -0.1 eV. The top and side views are shown in 
upper and lower panels, respectively. 
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